Nonlinear Collapse in the Semilinear Wave Equation in AdS 
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Previous studies of the semilinear wave equation in Minkowski space have shown a type of critical 
behavior in which large initial data collapse to singularity formation due to nonlinearities while 
small initial data does not. Numerical solutions in spherically symmetric Anti-de Sitter (AdS) are 
presented here which suggest that, in contrast, even small initial data collapse eventually. Such 
behavior appears analogous to the recent result of Ref. |T] that found that even weak, scalar initial 
data collapse gravitationally to black hole formation via a weakly turbulent instability. Furthermore, 
the imposition of a reflecting boundary condition in the bulk introduces a cut-off, below which initial 
data fails to collapse. This threshold appears to arise because of the dispersion introduced by the 
boundary condition. 
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Introduction: Recently, numerical studies of the gravi- 
tational collapse of a scalar field in asymptotically Anti- 
de Sitter (AdS) space found that the scalar field collapses 
to a black hole eventually for any initial amplitude for 
generic initial data [U [2] ■ This inevitability of black hole 
formation is understood as thermalization within the cor- 
responding conformal field theory (CFT) on the bound- 
ary of AdS according to AdS/CFT correspondence. The 
nature of the boundary of AdS is such that it can be 
reached in finite time and therefore the bulk is a bounded 
domain. As noted in Ref. pQ , there is evidence for similar 
nonlinear behavior in other, non-gravitational systems in 
bounded domains. 

It is in the context of such studies that one considers 
the dynamics of the semilinear wave equation in AdS 
for a complex scalar field <f> where p is an odd integer 
by adding a nonlinear potential term </> p to the scalar 
wave equation. Previous studies of this model (for a 
real scalar field) in Minkowski space have shown that 
for p = 7, the scalar field "collapses" to singularity 
formation for large initial data and disperses for small 
initial data [3J 0]. In AdS space, one might suspect 
that even small initial data, which would otherwise 
disperse in Minkowski space, would instead reflect off 
the AdS boundary repeatedly and eventually collapse. 
Indeed, numerical solutions presented here suggest that 
such behavior is present for the semilinear equation. 
Furthermore, placing a reflecting boundary at some 
finite radius places a lower limit on the size of collapsing 
initial data, similar to that found in the gravitating 
case [5]. Some speculation about the nature of this 
change in behavior is given. 

Implementation: Adopting essentially the same nota- 
tion and form as in [5], the metric is assumed to be that 
of spherically symmetric AdS 



ds 2 = 



-dt 2 + dx 2 



sin 2 x dQj 



(1) 



where I is the scale-size of the AdS spacetime, dVL 2 dl 
is the metric of S^ -1 . The domain extends from the 
origin, x = 0, to the boundary of AdS, x = n/2. In- 
troducing the auxiliary quantities $i = {8/ dx) (<pi) and 



ITj = (d/dt) ((f>i), one rescales according to 



d t <f>i 



n, 



Or 



cos u_± x 



(2) 
(3) 
(4) 

cos" _i x cos"~^ x 

Here, <j>i for i = 1,2 indicates the real and imaginary 
components of the scalar field. The second-order nonlin- 
ear wave equation becomes a system of three first-order 
equations 



= n 4 , 



sin d 1 x 



cos d ~ 1 xfl 



sin d 1 x - 
$ 



COS X 




(5) 
(6) 



(Hcos^ 1 



p-1 - 



The U(l) symmetry of the complex scalar field results in 
a conserved charge given by 

r /2 / - \ 

Q = J dx tan^xcos 2 ^ -1 ) x (ili^ - IIs^J . (8) 
The total energy E in the system is computed as 



E 



tt/2 



dx tan d_1 x cos 2(d ~ 1) x 



p+i 



p+i 

(9)' 



where sum over i is implied. The rescaled scalar quanti- 
ties are subject to the boundary conditions <j>i(ir/2,t) — 
0, $,(71-/2, <) = 0, and Ui(n/2,t) = 0. 

The initial data used in this work is generalized from 
that presented in [T] and is a subset of that used in [S] 

4>i(x,0) = —e-^^ cos^xS} (10) 



fk(x,0) = W>i(z,0), 



(11) 



where a, e, and u> are arbitrary constants. The initial 
profile for $ is determined by the spatial derivative of 
the initial scalar field 



&i(x, 0) = (1 - d) tan x qb t + cos 



(12) 
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FIG. 1: Demonstration of convergence through a bounce with 
e — 1 and w = 5. Runs with successively doubled resolutions 
are compared where the run indicated as "32" has a resolution 
smaller by a factor of 2 5 than the base resolution run "1." (top) 
The order of convergence for the field IIi as a function of 
time is shown. The other fields converge similarly. Runs were 
carried out with multiples as shown of a base resolution. The 
plot shows that the scheme is generally accurate roughly to 
fourth order, except when the pulse is at a boundary where it 
drops to third order, (middle) The fractional change in energy 
log 10 \E(t) — E(0)\ /E(0) as a function of time for the four 
highest-resolution runs. The conservation of energy improves 
rapidly with resolution and conserves energy to better than 
about one part in 10 s . (bottom) The fractional change in 
charge versus time. The conservation similarly improves with 
resolution. 



For w^O, the initial data is charged. 

As discussed in [3j 0] , when the amplitude of the scalar 
field becomes large, the nonlinear potential term domi- 
nates and focuses the scalar pulse, leading to singular- 
ity formation at some collapse time t c . Numerically, 
such collapse is "detected" by stopping the code when 
\4>\ > 10, where the value 10 is arbitrary. 

These equations are solved using the same infrastruc- 
ture and methods described in |T|, although no metric 
terms are required since the model evolves in pure AdS. 
Tests of the code indicate that it converges to better 
than third order in the grid spacing and conserves total 
charge and energy. The convergence order and fractional 
changes in energy and charge are presented for a typical 
run in Fig. [I] 

Results: Typical results are represented in Fig. [2] 
Charged (blue open circles), uncharged (cyan stars), and 
charged higher dimensional (magenta, solid squares) evo- 
lutions are shown for a wide range of e. The collapse 
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FIG. 2: Time to collapse as a function of e in the initial 
data. As e decreases, the initial pulse executes more bounces 
and the collapse time increases. Collapse occurs roughly 
in multiples of At — ir and collapse times roughly follow 
log 10 (t c ) = At + B. Compare to the bottom frame of Fig. 
1 of [2] which shows the time of collapse to black hole for a 
scalar field. Charged (blue, open circles) initial data follows 
a trendline with A = —1.11 and B — 5.91, while uncharged 
(cyan stars) configurations follow A — —1.08 and B = 5.78. 
A charged configuration in d = 4 (magenta solid squares) 
follows a slightly less negative trend with A = —0.914 and 
B = 5.14. 

times get progressively longer (roughly in multiples of 
the crossing-time n). The important point here is that 
the scalar executes an increasing number of reflections off 
the boundary with decreasing e and eventually collapses. 
As also found in [5] for the gravitating scalar, the global 
charge does not appear to have any significant effect. 

The collapse times for a different family of initial data 
are shown in Fig. [3] This family is characterized by a — 
1/8, twice the value characterizing the data shown in 
Fig. [2j These data similarly suggest scalar collapse is 
inevitable for any initial amplitude. 

Ref. [5] studied gravitational collapse on a restricted 
computational domain by imposing reflecting boundary 
conditions at some largest value x max of the coordinate 
x. In such a space, the scalar pulse does not "see" the 
full AdS background but nevertheless still evolves within 
a bounded domain. That the evolutions showed collapse 
after a number of reflections supported the view that the 
cause behind the inevitable collapse was that the domain 
was bounded, not that the space was AdS in particular. 
The boundedness of the domain allowed for the nonlin- 
earity of gravity to act continuously. However, Ref. [5] 
also reported that, although collapse did occur after a 
number of reflections, for small enough initial amplitude, 
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FIG. 3: Collapse times for families characterized by a = 0.125 
(in contrast to families with a = 0.0625 shown in Fig. [2j|. 
Initial data with rotation (blue circles) displays a transition 
from A = -1.55 and B = 6.38 to A = -1.96 and B = 7.87. 
Initial data without rotation (magenta crosses) follows the 
trendline described by A = —1.36 and B — 5.75. 
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FIG. 4: Collapse times for finite-domain evolutions. The same 
results for the ui = 5, d — 3 evolutions in the full domain 
from Fig. [2] are shown to facilitate comparison and contrast 
between full-domain and finite-domain results. Note that as 
the outer boundary is moved inwards, the minimum value of 
e that produces collapse decreases. For sufficiently small e, no 
collapse is detected. The collapsing solutions closest to the 
threshold have values of e within a few parts in a million of 
the e values of solutions which do not collapse. 



collapse did not occur, even after many reflections. 

In light of this interesting effect, evolutions of the 
nonlinear wave equation here are conducted within a 
restricted domain with the identical, reflecting bound- 
ary condition of Ref [S]. In particular, the conditions 
0i(^max ; t) = and IIi(x ma x, t) = are enforced. Studies 
with domains extending only to x max = 0.95,0.90,0.85, 
and 7r/4 w 0.785 are shown in Fig. |4j along with those 
within the full AdS spacetime (x max = 7r/2) for compar- 
ison. 

These restricted domain runs behave in stark con- 
trast with the full-domain results. Evolutions which col- 
lapse essentially immediately (before any reflection off 
the boundary) nearly coincide as would be expected. Af- 
ter one reflection, one would expect some range of e- 
values that would result in collapse, and therefore these 
restricted domains collapse faster than the full-domain 
evolutions simply because the geometric "concentration" 
at the origin happens earlier. However, as one con- 
siders even smaller e values, one sees that the collapse 
time grows apparently asymptotically to infinity. In 
other words, the evolutions suggest the existence of some 
threshold value of e* below which finite-domain evolu- 
tions no longer collapse. This minimum value appears to 
decrease with decreasing x max . 

As e is decreased in increasingly finer steps near the 
threshold, more bounces are observed before collapse, as 
indicated in Fig. [4] It is not clear yet what happens in the 
limit that e approaches some critical e* , and in particular 
whether the number of bounces continues to increase. 

That small-e, restricted-domain evolutions fail to col- 
lapse argues against the idea that the important effect 
of AdS is its introduction of reflections. Instead perhaps 
some other property of AdS accounts for the inevitable 
collapse and by restricting the domain, this property is 
affected. That is not the argument presented here. In- 
stead, there appears to be another effect introduced by 
the restriction of the domain, an effect not seen in the 
full domain case. 

Consider the wave equation with no nonlinear term, 
that is Eq. Q without the term. In spherically 

symmetric Minkowski space, the solutions take the sim- 
ple form r<p(r,t) — f(r ± t) where r is the radial coor- 
dinate and /() is some function. There is no dispersion 
of /(), although (f>(r,t) spreads out with increasing r. In 
AdS, however, things are not so simple and one expects 
that the metric terms will cause dispersion. In Fig. [5] 
the behavior of IIi(0,i) at the origin is shown at vari- 
ous times. In the full domain, the pulse "reflects" off the 
AdS boundary without inversion and flips when it im- 
plodes through the origin. Besides the inversion, the be- 
havior of the pulse remains the same, even at late times. 
The periodicity results from the fact that the modes of 
the scalar field in AdS have integer values, and hence a 
generic initial configuration will repeat with a period (at 
most) of At = 2ir 7 . In the restricted domains, the 
pulse inverts at both the outer boundary and the origin, 
and it disperses. Because the restricted domain does not 
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FIG. 5: Demonstration of dispersion with restricted domain 
evolutions of the "linear" wave equation. Shown is the be- 
havior of IIi at the origin versus time for evolutions in which 
the nonlinear term cf 7 has been removed. For reference is 
shown the full domain result (solid blue) which maintains its 
shape (except for an inversion of sign) with each reflection. 
The restricted domain evolutions show no inversion but do 
change shape. The bottom frame shows the late time behav- 
ior where the dispersion of the restricted domain evolutions is 
quite apparent. Also shown (dotted black) is the full domain 
pulse from t w 7r, shifted in time and inverted in sign. That 
it nearly overlays the late-time, full domain result indicates a 
lack of dispersion when evolving in the full AdS space. 



allow for integer eigenvalues, it is not periodic. 

The results suggest that this dispersion competes 



with the nonlinearity (be it gravitational or simply 
some scalar potential term). In this way, restricted 
domain evolutions of small initial data arc dominated 
by the dispersion whereas strong initial data are instead 
dominated by the nonlinearity. In the special case where 
the full domain is allowed, no dispersion occurs, and any 
initial data will eventually collapse. 

Conclusions: These studies provide an example of a 
non-gravitational, hyperbolic system which results in 
sharpening, and eventual collapse, of an initial pulse. 
This result appears analogous to the gravitational col- 
lapse of a scalar field in AdS. 

Evolutions within a restricted domain with reflecting 
boundary conditions also behave similar to the gravitat- 
ing scalar [5]. In particular, such evolutions collapse for 
only a limited range of decreasing amplitude and stop col- 
lapsing below some threshold initial size. The imposition 
of a reflecting boundary condition introduces dispersion 
where it would otherwise not appear, and this dispersion 
appears to compete with the weak turbulence pQ that 
transfers energy to higher frequencies. As one considers 
smaller initial data, one observes a transition from the 
dominance of weak turbulence to dispersion. 

It should be noted that not all scalar configurations 
in the gravitating case are unstable [6], and it would be 
interesting to study the p — 5 case of the semilinear equa- 
tion because, at least in Minkowski space, it possesses a 
static solution. 
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